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To determine the strength of essential electronic and magnetic interactions in the iridates S^IrCU 
and Ba2lr04 -potential platforms for high-temperature superconductivity- we use many-body tech- 
niques from wavefunction-based electronic-structure theory. Multiplet physics, spin-orbit interac- 
tions, and Ir-O hybridization are all treated on equal footing, fully ab initio. Our calculations 
put the lowest d-d excitations of Sr2lr04/Ba2lr04 at 0.69/0.64 eV, substantially lower than in 
isostructural cuprates. Charge-transfer excitations start at 3.0/1.9 eV and the magnetic nearest- 
neighbor exchange coupling is 51/58 meV. Available experimental results are fully consistent with 
these values, which strongly constrains the parametrization of effective iridate Hamiltonians. 
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Introduction — Electronic correlations in quasi two- 
dimensional (2D) materials, especially transition metal 
compounds, have been a central issue in condensed mat- 
ter physics since the discovery of the cuprate supercon- 
ductors. While electron-electron interactions are very 
substantial in 3d transition metal compounds [l|, they 
become progressively weaker when going to heavier tran- 
sition metal elements, i.e., Ad and 5d systems, as the d 
orbitals become more and more extended. Interestingly, 
the relativistic spin-orbit coupling (SOC) follows the op- 
posite trend - it increases progressively. In 5d transition 
metal compounds, e.g., iridates, the intriguing situation 
arises where these interactions meet on the same energy 
scale. The interplay between crystal-field effects, local 
multiplet physics, SOC's, and intersite hopping that sub- 
sequently emerges has opened up a new window of inter- 
est in strongly correlated compounds, offering novel types 
of correlated ground states (GS's) and excitations. 

The more exotic examples are possible topological 
states in iridates, such as a topological Mott insula- 
tor |2j, y] , a Weyl semimetal and axion insulator [4j , and 
the possible realization of the long-sought Kitaev model 
with bond-dependent spin-spin interactions [5Hl|- The 
2D square-lattice iridates such as Sr2ix04 and Ba2lr04 
are on the other hand appealing because of their per- 
ceived structural and magnetic similarity to La2Cu04, 
the mother compound of the cuprate high-T c (HTC) 
superconductors, which has promoted them to novel 
platforms on which HTC superconductivity may be de- 
signed [8h11|. To put such considerations on a firm 
footing it is essential to quantify the different coupling 
strengths and energy scales, as they for instance appear 
in effecti ve Ham iltonian descriptions of these correlated 
systems 
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Whereas the relevant single-particle hopping integrals 
can be extracted from tight-binding fits of density- 
functional-based band structures [8J, ll2Ml4J | . for a well- 
founded evaluation of interactions strengths and exci- 
tation energies from first principles one has to rely on 



wavefunction many-body approaches 15H17j|. Here, we 
present the outcome and discuss the implications of ab 
initio many-body calculations for the two so called 214 
iridates, S^IrCU and Ba2ir04. We find that for Sr/Ba 
214 the lowest d-d excitations of the magnetic GS dou- 
blet occur at 0.69/0.64 eV, substantially lower than in 
isostructural cuprates, where these excitations start at 
1.0-1.5 eV [lfjf. The charge-transfer (CT) continuum of 
the 2D iridates commences at 3.0/1.9 eV, with the low- 
est CT excitations confined to the out-of-plane oxygen 
ligands. It is quite remarkable that in spite of the very 
different Ir-O-Ir bond angles in the two materials, the 
strength of the magnetic nearest-neighbor (NN) exchange 
coupling is rather similar, 51 and 58 meV, respectively. 
This perhaps unexpected result is a direct consequence of 
the strong SOC and it sets the 2D iridates apart from the 
3d transition metal systems, in which superexchange cou- 
plings in general depend extremely sensitively on these 
geometrical factors. 

Computational method — For our calculations we rely 
on multiconfiguration self-consistent-field (MCSCF) and 




FIG. 1. Sketch of the 5-octahedra cluster used for the cal- 
culations of crystal-field and CT excitations in Sr2lr04. Ir, 
O, and Sr ions are shown in light blue, red, and pink, respec 
tively. In Ba2lrC>4, the Ir-O-Ir bonds are straight 
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multireference configuration-interaction (MRCI) tech- 
niques from modern quantum chemistry 20]. Multior- 
bital and multiplet physics, SOC's, and O 2p to Ir 5d 
CT effects are all treated on equal footing, fully ab ini- 
tio. Since spin-orbit interactions can be switched on and 
off by hand, we can filter out and explicitly determine 
the effect of SOC on the electronic structure and exci- 
tations. That the quantum-chemical calculations offer 
insights into the electronic structure of correlated solids 
that go substantially beyond standard density-functional 
approaches is well-established by now, e.g., for the 2D 
Cu oxides 15l-tl7|. The magnetic exchange constants 
[15j and d-d excitations [16[ obtained in those systems 
by quantum-chemical calculations are in excellent agree- 
ment with the experiment. Further, spin renormalization 
effects on the cuprate quasiparticle bands have also been 
determined 17|. Since the quantum-chemical method- 



TABLE I. Splittings within the Ir ti g shells in Sr2lr04 and 
Ba 2 Ir0 4 , with and without SOC. MCSCF and MRCI results 
for 5-octahedra embedded clusters, hole representation. 



ology does not use any ad hoc parameters, such calcu- 
lations are very different from, e.g., hybrid schemes jl 31 ] 
based on density-functional theory (DFT) and dynamical 
mean-field theory (DMFT) and provide information on 
the electronic structure and excitations that is comple- 
mentary to for instance a DFT+DMFT approach [13J . 

In the spirit of modern multiscale electronic-structure 
approaches, a given region of the extended crystal is 
treated by advanced quantum-chemical many-body tech- 
niques while the remaining par t of the solid is modeled 
at a more elementary level 15l4l7|. For each of the two 
compounds, we first designed 5-octahcdra clusters (see 
Fig. 1) to address the nature of the GS electron configu- 
ration and of the lowest crystal-field (CF) and CT excited 
states. The ten Sr or Ba ions next to the central Ir06 oc- 
tahedron are also included in these clusters. The remain- 
ing part of the crystal is modeled as a large array of point 
charges (PC's) which reproduce the Madelung potential 
in the cluster region. For the GS calculations, the orbitals 
within each finite cluster are variationally optimized at 
the MCSCF level. All Ir t2 g functions are included in the 
active orbital space [20[, i.e., all possible electron occu- 
pations are allowed within the ti g set of orbitals. On-site 
and CT excitations are afterwards computed just for the 
central IrO@ octahedron while the occupation of the NN 
Ir valence shells is held frozen as in the GS configura- 
tion. We analyze then the correlated electronic structure 
of these compounds at the MCSCF and MRCI levels of 
theory, both with and without SOC's [UHl]. To extract 
the NN superexchange J, we employed PC-embedded 8- 
octahedra clusters. For those clusters, high-spin and low- 
spin magnetic configurations are computed for two active 
Ir sites. To reduce the complexity of the problem and the 
computational effort, the six adjacent Ir 4+ ions are mod- 
eled as closed-shell Pt 4+ t% g species. 

All calculations were performed with the molpro 
quantum-chemical software [23j. We used energy- 
consistent relativistic pseudopotentials for Ir, Sr, and 
Ba 2J, |25[ and Gaussian-type basis functions from the 
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MOLPRO library. Basis sets of quadruple-zeta quality [20j 
were applied for the valence shells of the central Ir ion 
and triple-zeta basis sets for the ligands of the central 
octahedron and the NN Ir sites. For the central Ir ions 
we also used two polarization / functions. For farther 
ligands around the NN Ir sites we applied minimal basis 
sets. All occupied shells of the Sr 2+ and Ba 2+ ions were 
incorporated in the large-core pseudopotentials 25[ and 
each of the Sr 5s and Ba 6s orbitals was described by a 
single contracted Gaussian function. The ligand bridg- 
ing the two active Ir sites in the 8-octahedra clusters was 
described with quintuple-zeta valence basis sets and four 
polarization d functions. 

Splittings of Ir t2 g levels — For the (formal) Ir + t|n 
valence electron configuration, the strong SOC tends to 
split off a low energy j = 1/2 doublet. Much attention 
has been recently paid to obtaining the size of this split- 
ting, the most basic information on the SOC strength 
A [8j, Il4[ |26 1 , and the nature and size of the intersite spin- 
exchange interactions in a number of different iridates 



[a-l7l lKx l27H29t ]. As concerns the magnetic interactions, 
it has been pointed out that fine details such as the size 
of the splitting(s) within the valence Ir t2 g shell can be 
extremely relevant. In the layered honeycomb systems, 
for instance, the trigonal CF splittings ap pea r to strongly 
influence the nature of the magnetic GS [28j, [29( . 

Whereas the crystal structures of the 214 iridates are 
similar to that of La2Cu04, an important difference be- 
tween the Sr and Ba 214 systems is that in the former 
material the IrOg octahedra are rotated about the c axis 
such that the in-plane Ir-O-Ir angles are not 180° while 
in the latter the Ir-O-Ir bonds are straight. In both 
Sr2lr04 and Ba2ir04, the apical Ir-0 bonds are longer 
than the in-plane bonds, 2.06 vs. 1.98 A in Sr2lr04 and 
2.15 vs. 2.01 A in Ba 2 Ir0 4 [lElll- Results for the split- 



tings within the Ir ti g shells in the two materials are 



TABLE II. Nearest-neighbor exchange coupling J in Si^IrCU 
and Ba2lr04. Negative/positive values denote FM/AF in- 
teractions. MCSCF and MRCI results for clusters with two 
active magnetic centers, see text. 
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listed in Table I. Interestingly, the calculations in which 
the SOC is switched off show that, in the electron repre- 
sentation, the lowest, doubly occupied ti g levels are d xz 
and d yz in Sr2lr04 and d xy in Ba2lr04. By MRCI cal- 
culations, the splittings within the ti g shell are 0.06 and 
0.27 eV, respectively (see Table I). We included in the 
MRCI treatment single and double excitations from the 
Ir tig and all O 2p orbitals at the central octahedron. 

Test calculations in which for the crystal structure of 
Sr2lr04 we replaced the Sr ions by Ba show that the 
d xz and d yz levels are still the lowest, as for the set of 
results listed in the upper part of Table I. Also, if we use 
the crystal structure of Ba2lrC-4 but replace Ba by Sr 
we still find that the lowest-energy ti g level is d xy , as 
for the results given in the lower part of the table. This 
shows that our findings of a different energy order of the 
Ir tig levels for Sr and Ba 214 in absence of SOC are 
strictly related to local Ir-0 bonding. Nevertheless, when 
SOC is accounted for, the splittings between the local 
doublet j = l/2 and quartet 7=3/2 states are about the 
same in the two compounds. As mentioned above, these 
calculations were performed on embedded 5-octahcdra 
clusters. In a nonrelativistic picture, the five Ir S = 1/2 
sites can couple to sextet, quartet, and doublet states. 
All those different spin configurations entered the spin- 
orbit calculations so that the energy windows of 0.65- 
0.81 cV in Sr 2 Ir0 4 and 0.69-0.82 eV in Ba 2 Ir0 4 act ually 
include a total number of 64 spin-orbit final states [30j. 
A simplified picture can be obtained by replacing the 



Ir 4 + tL NN's by closed-shell Pt 4 + t 6 
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2g ions. The energies 



of the two components of the j = 3/2 quartet, split off 
by the tetragonal field, are in that case 0.69 and 0.73 in 
Si'2lr04 and 0.64 and 0.71 in Ba2lr04. These results for 
the energy separation between the j = 1/2 and j = 3/2 
states imply an effective spin-orbit coupling parameter A 
for Sr/Ba 214 of 0.46/0.43 eV, which is surprisingly close 
to values of 0.39-0.50 eV extracted for Ir 4+ impurities 
from electron spin resonance measurements 3l|, |32j . 

CF and CT excitations — MRCI calculations without 
SOC put the CF excitations from the tig to the e g levels, 
t\ g ^ft\ g e\ and t\ g ^tt\ g e\, between 2.8-6 eV in Sr 2 Ir0 4 
and between 2.0-5.0 eV in Ba2lr04. The lowest O 2p to 
Ir 5d CT states in MRCI calculations without SOC are 
at 3.6 eV above the GS in Sr2lr04 and at 1.7 eV with 



respect to the GS in Ba 2 Ir04. When SOC is switched 
on, the CF and CT configurations mix and the lowest 
of those spin-orbit states require excitation energies of 
3.0 eV in Sr 2 Ir0 4 and 1.9 cV in Ba 2 Ir0 4 . For Sr 2 Ir0 4 , 
those values fit well the resonant inelastic x-ray scatter- 
ing (RIXS) [33| and optical absorption data 34j, with 
both types of measurements displaying marked spectral 
features between 2.5 and 8 eV. 

For both Sr and Ba 214, the computed CT gaps are 
relatively small. Nevertheless, the interaction and mix- 
ing of non-CT t^g and CT t% g L_ configurations via SOC 
is not substantial. It should be noted, however, that a 
major difference with respect to the 2D 5=1/2 cuprates 
is that in the iridates the lowest CT states have apical 
hole character, whereas in the 2D Cu oxides the lowest 
electron-removal states have predominant in-plane O 2p 
character, see, e.g., Ref. [171 ]. 

Magnetic exchange couplings — Our MRCI calcula- 
tions for the superexchange coupling J include single and 
double excitations from the Ir ti g and bridging-ligand 
2p orbitals. Quantum chemical calculations of magnetic 
couplings in terms of restricted configuration-interaction 
expansions involving the d orbitals at the active magnetic 
sites and p functions of the bridging ligand were previ- 
ously reported, e.g., for La 2 Cu04 35], providing J's that 
are in good agreement with experiment. All possible oc- 
cupations were allowed within the set of ti g orbitals at 
the two Ir sites in the present MCSCF calculations, which 
gives rise to nine singlet and nine triplet states in absence 
of SOC. These MCSCF wavefunctions were expressed in 
terms of orbitals optimized for an average of singlet and 
triplet states. All those eighteen states entered the spin- 
orbit calculations, both at the MCSCF and MRCI levels. 
J then corresponds to the difference between the ener- 
gies of the lowest spin-orbit singlet and triplet state. The 
three components of the lowest jtotai = 1 triplet display in 
both compounds rather small splittings of not more than 
3 meV, which shows that the mapping of our ab initio 
results on a symmetric Heisenberg model is controlled. 
The small effects related to the presence of anisotropic 
non-Heisenberg terms in the interactions will be the sub- 
ject of future work. 

The resulting J's, as shown in Table II, have a clear 
physical interpretation. Since the splittings within the 
tig shell are relatively small, the three different t\ con- 
figurations (see Table I) appear with similar weights in 
the spin-orbit GS wavefunctions. Careful analysis of the 
wavefunctions for clusters with two active Ir sites however 
shows that the different order of the ti g levels in Sr and 
Ba 214 gives rise to subtle differences with respect to the 
superexchange interactions. For corner-sharing plaque- 
ttes or octahedra and ti g active orbitals, the Anderson 
superexchange implies d-p-d 7r-type links. For Sr2lr04, 
when the SOC is switched off, the relevant hole orbital is 
the in-plane d xy component, as shown in Table I. Due to 
rotation of the IrOg octahedra and bending of the d xy - 



Px/y-dxy bond in Sr 214 [18[ the antiferromagnetic (AF) 
superexchange is not effective and the spin interaction J 
is actually ferromagnetic (FM), see Table II. With SOC, 
the tig hole acquires d xz and d yz character too. Because 



the d 



xz/yz Pz U'xz/yz 



orbital overlap is rather unaffected 



by bending of the Ir-O-Ir bond, the d-p-d superexchange 
interactions are much more effective in the presence of 
SOC. Consequently the magnetic exchange J turns AF, 
both in MCSCF and MRCI (2nd column in Table II). 

Ba2lr04, on the other hand, lacks a bending of the Ir- 
O-Ir bond [19j. The nature of the hole orbitals in the 
GS wavefunction in absence of SOC therefore matters 
less in this case: J is large and AF, also in absence of 
SOC's (3rd column in Table II). At the MRCI level and 
including SOC's, J turns out to be somewhat larger in 
Ba 2 Ir04 as compared to Sr 2 Ir04, 58.0 vs. 51.3 meV. 

Conclusions — Our detailed investigation of the elec- 
tronic structure of the 2D j = l/2 Ir oxides Sr2lr04 and 
Ba2ir04, with ah initio many-body quantum-chemical 
calculations, shows unequivocally that the peak feature 
observed at about 0.6 eV in the RIXS HSHHj on Sr 214 is 
directly related to transitions between Ir t| g j = l/2 and 
j = 3/2 states that are in an energy window of 0.65-0.81 
cV in Sr 214 and 0.69-0.82 eV in Ba 214. It is interesting 
to note that the peak at this energy in optical absorp- 
tion measurements was instead assigned to an excitation 
between Hubbard bands [34|. The experimentally ob- 
served "gap" between the 0.6 eV feature and higher ex- 
cited states that appear above 2 eV in Sr 214 [3, [33J, |34| is 



fully consistent with our calculations. We find that the Sr 
and Ba 214 compounds differ significantly with respect 
to their tetragonal crystal-field splittings. While in Sr 
214 the lowest t-x g states have xz and yz symmetry, in 
Ba 214 the lowest ti g levels are the xy orbitals. The cor- 
responding tetragonal splittings are 0.06 eV in Sr 2 Ir04 
and 0.27 eV in Ba2lr04. The calculations further show 
that without SOC the distorted Ir-O-Ir bonds in Sr 214 
would actually give rise to a FM nearest-neighbor mag- 
netic exchange J. However, the strong SOC reverses the 
sign of J and causes antiferromagnetism. We find that 
the dominant magnetic exchange interaction is of Heisen- 
berg type and for Sr 214 has a substantial coupling con- 
stant J = 51 meV, which compares well with the value of 
45 meV estimated from experimental data [5j and with 
RIXS measurements that indicate J = 60 meV for this 
system [lOj. For Ba 214, we find a J that is even some- 
what larger, which renders it roughly a factor of 2 smaller 
than the J's in typical 2D cuprates [15j. This might in 
itself be encouraging for a scenario of magnetic medi- 
ated superconductivity in doped 214 iridates. However, 
the fact the lowest charge-transfer states in 214's have 
instead of planar, apical oxygen character, points at a 
destabilization of Zhang-Rice singlets and consequently 
an effective low energy Hamiltonian that is different from 
the standard t-J model of doped cuprates. 
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